High-Energy Neutrino Signatures of Dark Matter Decaying into Leptons 
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Decaying dark matter has previously been proposed as a possible explanation for the excess high 
energy cosmic ray electrons and positrons seen by PAMELA and the Fermi Gamma- Ray Space Tele- 
scope (FGST). To accommodate these signals however, the decays must be predominantly leptonic, 
to muons or taus, and therefore produce neutrinos, potentially detectable with the IceCube neutrino 
observatory. We find that, with five years of data, IceCube (supplemented by DeepCore) will be 
able to significantly constrain the relevant parameter space of decaying dark matter, and may even 
be capable of discovering dark matter decaying in the halo of the Milky Way. 
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A number of recent measurements of cosmic ray elec- 
trons and positrons have been interpreted as possible evi- 
dence for dark matter |H, 0, 0, 0, 0, E| . In particular, it has 
been suggested that the features of the e + e~ spectrum 
and positron fraction reported by PAMELA [1], ATIC 
0, PPB-BETS0, the Fermi Gamma Ray Space Tele- 
scope (FGST) [JHEAT @, and AMS-Oljl may origi- 
nate from either DM annihilations 0, 0, 0, [l0, [U, E2] or 
decays Q taking place in the Galactic Halo. It is also 
possible that these signals originate from astrophysical 
sources, rather than from new high energy physics [151 ]. 

Although annihilating dark matter could potentially 
generate the observed anomalous cosmic ray features, at- 
tempts to do so face a number of challenges. Firstly, 
the spectrum of electrons and positrons predicted to be 
generated in the annihilations of most dark matter candi- 
dates is much too soft to fit the observations of PAMELA 
and FGST [1, 0|. If WIMPs annihilating throughout 
the halo of the Milky Way are to produce the spectral 
shape observed by these experiments they must annihi- 
late mostly to charged leptons. While models have been 
proposed in which this is the case [Io|, [H| , many of the 
most often studied WIMP candidates (including MSSM 
neutralinos) are predicted to annihilate dominantly to 
quarks and/or gauge bosons [16[. Furthermore, annihi- 
lations to non-leptonic final states tend to produce more 
cosmic ray anti-protons than are observed [17[ . Secondly, 
the dark matter annihilation rate that is required to gen- 
erate the observed spectrum of cosmic ray electrons and 
positrons is considerably higher than is predicted for a 
typical thermal relic distributed smoothly throughout the 
Galactic halo. To normalize the annihilation rate to the 
PAMELA and FGST signals, we must require either large 
inhomogeneities in the dark matter distribution which 



lead to a considerably enhanced annihilation rate (i.e. a 
"boost factor"), and/or dark matter particles which pos- 
sess a considerably larger annihilation cross section than 
is required of a thermal relic. This latter requires either 
a non-thermal production mechanism in the early uni- 
verse, or an enhancement of the annihilation cross sec- 
tion at low velocities, such as through the Sommerfeld 
effect pH EH or Breit-Wi gner enhancement [13|- In this 
light of these challenges, the observations from PAMELA 
and FGST are extremely surprising and pose an interest- 
ing challenge to the usual WIMP paradigm. 

In this paper we consider decaying dark matter as a 
possible origin of these excess electrons and positrons, 
as has previously been discussed in Refs. [3, Hfll- 
The only properties of the dark matter relevant to this 
work are its mass, its lifetime, and its decay channels. 
Similar to the discussion in the previous paragraph, any 
explanation of the PAMELA and FGST data requires 
preferentially leptonic decay products. Beyond these few 
phenomenological properties, the nature of the decaying 
particles is not relevant to our study. 

In order to confirm a particle physics origin of the 
PAMELA and FGST data, one would hope to observe 
other final state particles in addition to electrons and 
positrons. It has been shown that dark matter with a 
mass ~ 10 2-5 GeV that annihilates/decays into 
r~r + , or can reproduce the observed cosmic 

ray features @. In such a scenario, we also expect the 
associated production of gamma rays and neutrinos from 
the decay of the heavy charged leptons. 

In a previous letter [21], the authors examined the sen- 
sitivity of the IceCube detector to neutrinos from dark 
matter pair-annihilation in the inner Milky Way. As 
the annihilation rate rises with the dark matter density 
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squared, the inner galaxy is a promising region to observe 
such signatures. Unfortunately, at the Antarctic location 
of IceCube, the Galactic center is overhead, causing the 
signal to be swamped with background from atmospheric 
muons. However, the Deep Core extension of IceCube (to 
be completed in 2010) [22|, HH, the heavily instrumented 
inner volume of the detector can make use of the re- 
maining IceCube volume to reject the muon background, 
allowing for neutrino-induced showers to be identified. 
This reduces the background to those events from atmo- 
spheric neutrinos, which is not overwhelming compared 
to the expected signal. Using this technique, it will be 
possible to place bounds on the dark matter annihilation 
rate that are comparable to those required to explain the 
PAMELA/FGST results. 

A natural extension of this idea is to consider dark 
matter decay. In this letter, we consider dark matter 
decaying into muon pairs, tau pairs, or four muon final 
states, each of which can adequately explain the observed 
electron/positron spectrum, for dark matter lifetimes on 
the order of r - lO 26 " 27 seconds 0, H, Hj| . In all three 
of these cases, the subsequent decays of the muons and 
taus will create copious numbers of high energy neutri- 
nos, potentially detectable at IceCube. 

The flux of neutrinos from dark matter decay in the 
inner Milky Way is given by 
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Here, r = r is the decay width of the dark matter, m x 
is the dark matter mass, and dNi/dE is the differential 
flux of neutrinos of flavor i resulting from the decay. The 
dark matter distribution integrated over the line-of-sight 
over a solid angle Aft is given by 
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p x (s)ds ; J(An) = —- PSF*JdQ 
AS 2 J A q 

(2) 

where PSF is the point spread function of the instrument 
and p x (s) is the dark matter mass density distribution. 
Note that this differs from the definition of J used in 
Ref. [21], in which J oc p 2 . As we are interested here 
in decay rather than annihilation, only one power of the 
density appears. 

The spectrum of neutrinos from the decays depends 
on the mass of the dark matter and on the dominant 
decay channel. For muon-channel decays (p + p~ and 
p + p~ p + p~), the muons themselves decay to v^ev e . For 
decays through taus, there are many final states avail- 
able, including r — > /i^ M ^ r , ev e v Tl as well as from the 
hadronic decays r — > 7r^ r , Kv T , ttttv t ^ and ttttttv t {24[. 
For this work, the spectrum of neutrinos from tau decays 
as a function of m x was determined numerically using 
Pythia [IK . It should be noted that for the four muon fi- 
nal state, it is assumed that the two p + p~ pairs originate 
from the decays of back-to-back parent particles. As a 



result, the rest energy m x is split evenly between the two 
pairs, rather than distributed among the four particles as 
per naive phase- space. 

The primary backgrounds for this observation consist 
of atmospheric muons and neutrinos. As addressed pre- 
viously, the IceCube detector itself can be used to veto 
muons inside of the volume of DeepCore, leaving only 
neutrino-induced showers to compete with. For the spec- 
trum of atmospheric muon neutrinos, we use the results 
of Ref. [26[ , which are in good agreement with the mea- 
surements of AMANDA [23. 1 

The effective area of the detector for neutrinos can be 
defined as 



A(E)np ice N A a vN (E)V(E), 



(3) 



where p ice = 0.9 g/cm 3 , N A = 6.022 x 10 23 g" 1 (to con- 
vert grams to nucleons), ct u n(E) is the neutrino- nucleon 
cross-section [28| and V(E) « 0.04 km 3 is the effective 
volume of the DeepCore detector for a neutrino-induced 
shower of energy E [231 ]. 

The directional capability of IceCube/DeepCore for a 
neutrino-induced shower above 1 TeV is expected to be 
on the order of 50° [23[. We conservatively consider the 
signal and background over a solid angle corresponding 
to a full half of the sky (2tt sr), acknowledging that our 
results would be strengthened if better angular resolu- 
tion could be obtained. Using an NFW profile, we inte- 
grate the dark matter distribution in the direction of the 
Galactic center over this solid angle. We take the energy 
resolution of the detector to be log(E max /E m i n ) ~ 0.3 

& 

For dark matter masses between 100 GeV and 30 TeV, 
we calculate the lifetime for which IceCube/DeepCore 
would provide either a 2<j exclusion limit or a 5cr discov- 
ery after 5 years of observation. The limits for the three 
decay channels under consideration are shown in Table H 
In Fig. [U we overlay our results on top of the regions of 
m x vs. r parameter space preferred by PAMELA/FGST. 
Also shown are the regions of parameter space excluded 
by HESS and VERITAS 0, 0, [p, 0, O as well as 
the limits from Super-Kamiokande [36|, [37[, FGST obser- 
vations of gamma rays from inverse Compton scattering 
below 10 GeV, and radio observations of the Galactic 
Center. In the figure, our results (blue and red 2 and 5cr 
lines from IceCube/DeepCore) have been superimposed 
upon a plot taken from Ref. [19[ with permission of the 
authors. In this plot, the dwarf spheroidal (dS-7) limit 
was originally determined in Ref. [341 ]. 



1 As the background is considerably larger than that from ^ e 's, 
any discrimination between electromagnetic and hadronic show- 
ers could be used to further reduce the backgrounds and improve 
the statistical reach of IceCube/Deepcore to the signal described 
in this paper. 
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FIG. 1: Exclusion limits in the dark matter lifetime versus mass plane, for decays to 



and . The green 



region is preferred by the observations of PAMELA, light red is preferred by FGST, and the dark red region is preferred by 
the combination of these measurements. Here, we assume that the excess of the cosmic e + + e~ spectrum measured by FGST 
beyond the simple power law is the dark matter decay signal (For alternative interpretations, see, e.g., [35]). The blue line 
labeled "2<j" is the projected 2a limit for 5-years of running at IceCube/DeepCore, while the red "5cr" line is the 5a discovery 
reach, assuming an NFW profile (the results change only mildly if another profile is adopted). The black line labeled "Super-K" 
is the limit from neutrino observations at Super-Kamiokande l36l |37ll. HESS bounds from the Galactic Ridge ("GR-7") [3ll ], 
Galactic Center ("GC-7") [30l ], and dwarf spheroidals ("dS-7") 32, 33] are also displayed. The dS-7 limit originally calculated 
in [3j]. Inverse Compton scattering limits from FGST are labeled as "IC," and radio observations of the Galactic Center are 
"GC-radio" 38]. Our results (blue 2a and red 5cr lines from IceCube/DeepCore) have been superimposed upon a plot taken 
from Ref. [19|| with permission of the authors. 



A noticeable dip in sensitivity is found for the case of a 
~100 GeV dark matter particle decaying to four muons. 
In this case, the maximum energy carried away by each 
muon is not much larger than the energy threshold of 
IceCube/DeepCore, approximately 20 GeV. As a result, 
most of the neutrinos from such decays are unobservable. 

From the results of Fig. HJ we see that IceCube, in- 
cluding the DeepCore extension, can probe the most 
interesting regions of the m x -r parameter space with 
five years of data after the planned completion in 2010. 
Due to its vastly larger volume, IceCube can probe 
dark matter decay lifetimes that are orders of magni- 
tude greater than those excluded by Super-Kamiokande. 
IceCube/DeepCore will be competitive with or more con- 
straining than inverse Compton and 7-ray measurements, 
and will serve as a complementary test of the dark matter 
interpretation of PAMELA and FGST results. 
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m x 
(GeV) 


Bin Size 
(GeV) 


5cr Detection 
(x -> M + M~) 


2cr Exclusion 

(x -> M + M~) 


100 


20 - 250 


r < 0.11 x 10 26 s 


r < 0.28 x 10 26 s 


150 


20 - 250 


r < 0.15 x 10 26 s 


r < 0.38 x 10 26 s 


300 


20 - 250 


r < 0.18 x 10 26 s 


r < 0.45 x 10 26 s 


500 


100 - 500 


r < 0.26 x 10 26 s 


r < 0.65 x 10 26 s 


1000 


150 - 800 


r < 0.40 x 10 26 s 


r < 1.0 x 10 26 s 


2000 


300 - 1500 


r < 0.55 x 10 26 s 


r < 1.4 x 10 26 s 


10000 


1500 - 8000 


r < 1.0 x 10 26 s 


r < 2.7 x 10 26 s 


30000 


5000 - 25000 


r < 1.4 x 10 26 s 


r < 3.9 x 10 26 s 


m x 
(GeV) 


Bin Size 
(GeV) 


5cr Detection 
(X -> r + r~) 


2a Exclusion 

(X-T+T-) 


100 


20 - 250 


r < 0.099 x 10 26 s 


r < 0.25 x 10 26 s 


150 


20 - 250 


r < 0.13 x 10 26 s 


r < 0.34 x 10 26 s 


300 


20 - 250 


r < 0.17 x 10 26 s 


r < 0.43 x 10 26 s 


500 


100 - 500 


r < 0.23 x 10 26 s 


r < 0.57 x 10 26 s 


1000 


150 - 800 


r < 0.34 x 10 26 s 


r < 0.87 x 10 26 s 


2000 


300 - 1500 


r < 0.47 x 10 26 s 


r < 1.2 x 10 26 s 


10000 


1500 - 8000 


r < 1.0 x 10 26 s 


r < 2.7 x 10 26 s 


30000 


5000 - 25000 


r < 1.2 x 10 26 s 


r < 3.3 x 10 26 s 


m x 
(GeV) 


Bin Size 
(GeV) 


5cr Detection 

(x -> mV"mV") 


2cr Exclusion 

(x -> mV~mV~) 


100 


20 - 250 


t < 0.0081 x 10 26 s 


r < 0.021 x 10 26 s 


150 


20 - 250 


r < 0.068 x 10 26 s 


r < 0.17 x 10 26 s 


300 


20 - 250 


r < 0.11 x 10 26 s 


r < 0.27 x 10 26 s 


500 


20 - 250 


r < 0.18 x 10 26 s 


r < 0.44 x 10 26 s 


1000 


100 - 500 


r < 0.25 x 10 26 s 


r < 0.65 x 10 26 s 


2000 


150 - 800 


r < 0.40 x 10 26 s 


r < 1.0 x 10 26 s 


10000 


600 - 3000 


r < 0.82 x 10 26 s 


r < 2.1 x 10 26 s 


30000 


1500 - 8000 


r < 1.2 x 10 26 s 


r < 3.2 x 10 26 s 



TABLE I: Limits on lifetime, r, of dark matter decaying into 
/x + /i - (top), t + t - (middle), or /x + /i - /i + /x - (bottom) after 5 
years of data collection at IceCube. Limits are placed for a 
discovery at 5cr as well as a bound for 2a exclusion. 
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